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Q I We present results from 1078 high resolution spectra of 990 stars in the young 

^ I open cluster NGC 2264, obtained with the Hectochelle multiobject echelle spec- 

ie I trograph on the 6.5m MMT. We confirm 471 stars as members, based on their 

radial velocity and/or Ha emission. The radial velocity distribution of cluster 

O 



members is non-Gaussian with a dispersion of a ~ 3.5km s . We find a substan- 
tial north-south velocity gradient and spatially coherent structure in the radial 



^ I velocity distribution, similar to that seen in the molecular gas in the region. Our 

j^ I results suggest that there are at least three distinguishable subclusters in NGC 

2264, correlated with similar structure seen in ^^CO emission, which is likely to 

CS I be a remnant of initial structure in this very young cluster. We propose that this 

Q ' substructure is the result of gravitational amplification of initial inhomogeneities 
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cn 



during overall collapse to a filamenta r y dist ribution of gas and stars, as found in 



O ' simulations by iBurkert &: HartmannI (120041 ) 
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Introduction 



Studies of young galactic clusters associated with molecular clouds can provide insight 
into the early history and formation of these objects. One of the most famous and best 
studied young cluster is NGC 2264. Many investigations of this region have been undertaken, 
building on the early work of Herbig (1954). A number of studies suggest low extinction 
toward NGC 2264 {E{B -V) = 0.082 by Walker (1956); E{B -V) = 0.071 ±0.033 by Sung 
et al. (1997); E{B -V) = 0.066 ±0.034 by Park et al. (2000)), although the large amount of 
molecular gas in the region (§2.4) suggests that a substantial number of sources embedded 
in dust remain to be found. The optical nebular emission is less intense than in the Orion 
Nebula Cluster (ONC), so contamination of stellar spectra from nebular emission is less of a 
problem in NGC 2264. There is a dark cloud immediately behind the cluster which obscures 
most of the b ackground stars, thus reducing the number of potential impostor members 
jHerbieJEgsi ). Thus, despite the fact that NGC 2264 is more distant than the ONC and 
other nearby regions of star formation ( 800 pc by Walker 1956; 800 pc by Sagar & Joshi 
1983; 950 ± 75 pc by Perez et al. 1987; 760 ± 40 pc by Park et al. 2000) the low extinction 
and weak nebular emission make NGC 2264 a favorable and important target for star and 
cluster formation studies. 

Age estimates of NGC 2264 have been reduced as the qu ality of NGC 2264 observations 
have improved. Ear ly estin iates suggested a ges 3 -5 Myr (jHarrid 1 19761 : iPark et al.l 12000 



Palla &: Stahler 



Ramirez et al. 



19991) . while iFlaccomio et al.l (2000) estimate an age closer to 1-2 Myr, and 
(I2OO4J ) found that most of the X-ray stars were younger than 3 Myr. A 



major reason for the discrepancy betwe en different estima tes is the differing weight placed on 
specific mass ranges. As pointed out by iHartmanru (120031 ). ages inferred fr om 1 — 2 Mq stars 



in young associations tend to be older than those of the lower-mass stars. IHartmanru (|2003l ) 
argued that the ages of these slightly higher-mass stars are systematically overestimated by 
neglecting "birthline" effects (e.g., pre-r nain sequence contractio n fror n finite starting radii); 
this leads us to favor the lower ages of IFlaccomio et al.l (120001 ). and iRamirez et al.l (120041 ) 
based on the lower mass stars, which are more consistent with the large amount of molecular 
gas associated with the cluster. 

The spatial distribution of the stars in NGC 2264 is extended alo ng the north-sout h 
direction, paralelling the extension of the molecular gas in the region (JRidge et al.ll2003l ). 
This feature, consistent with the general youth of the cluster, suggests that NGC 2264 may 
still exhibit the morphology of the gas cloud from which the protostars coalesced. The 
early evolution of stellar systems may be strongly affected by the density and dynam ics of 
the environments in which they form. Numerical simulations (IScally fc Clarkd 120021 ) show 
that the lifetime of any primordial substructure is short, a few million years or less. Radial 
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velocity studies may be able to provide clues to cluster formation for very young regions by 
revealing coherent patterns in the motion of sub-regions. 

In this paper we report the initial results of our continuing radial velocity (RV) survey 
of NGC 2264, in the search for kinematic signs of primordial structures in young stellar 
clusters. In §2 we describe the observations and data reduction; in §3 we compare the radial 
velocity distribution of the stellar and gaseous component, the latter as traced by ^^CO 
measurements. Finally, we discuss possible implications of our results for cluster formation 
scenarios. 



2. Observations and data reduction 
2.1. Target selection 

We drew our target sample from several sources, starting with surveys of Ha emission- 



l ine objects (Her big 1954; Reipurt 



fISung et al.lll997l . l2004l : IPark et al. 



letal. 2004 - R04 hereafter), narrow-band photometry 



20001). and samples of youn g objects selected based upon 



their IR and U V excess (JLada et al.l Il993l : iRebuU et al.l |2002| ) . All of these surveys have a 
bias toward finding (accreting) Classical T Tauri stars (CTTS) and against finding pre-main 
sequence (PMS) stars with weak or no disks. Therefore, by supplementing the se samples with 
a large catalog of X-ray sources (JFlaccomio et al.ll2000l ; iRamirez et al.ll2004j ). also sensitive 
to weak-emission (non-accreting) T Tauri stars (WTTS), the bias towards accreting stars 
can be corr ected. These surve ys, together with proper motion (jVasilevskis et al.lll965l ) and 
variabihty (JLamm et al.ll2004l ) surveys provide a very good input catalog of potential cluster 
members. 

The se lected targets c o me fr om three of the catalogs above: a substantial fraction, 1389 



stars, from iRamirez et al. 



stars from Park et al 



j 20041 ) X-ray sources; 408 Ha emission stars from R04; and 100 



(I2OOOI ). As some of the stars appeared in more than one of these three 



cataloges, the final target list contained 1563 stars. 



2.2. Spectroscopic observations 



The stellar spectra were obtained with Hectochelle (ISzentgyorgyi et al.lll998l ). a fiber- 
fed, bench mounted echelle spectrograph that operates at the post-conversion MMT. The 
MMT is a 6.5m Cassegrain teles cope and the Hectoch elle is designed to operate with the 
MMT in its f/5, wide field mode (JFabricant et al.ll2004l ). The robot positioner places all 300 
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fibers buttons for observation in five minutes with an accuracy of 25yum or better. The fiber 
diameter is 250/im, or 1.5" at the MMT f/5 plate scale. The resolution of the Hectochelle 
is ~34,000 and the the peak efficiency is 7% of the light striking the primary mirror. Be- 
cause the Hectochelle and Hectospec, a moderate dispersion fiber fed spectrograph, both 
use the same fiber feed, but have different magnifications, only 240 fiber images fall on the 
Hectochelle focal plane format, while Hectospec exploits all 300 available fibers. 

The preferred order for radial velocities with Hectochelle is centered at 5225A and 
includes the Mg triplet. Unfortunately, a failure of the original coatings (now replaced) on 
the camera and collimator mirrors of Hectochelle reduced the refiectivity in this spectral 
region by a factor of up to 6 per mirror. Also, our targets have brighter continuua at Ha 
than at 5225A, helping to compensate for fewer lines. Therefore we decided to use the 190A 
wide order centered at Ha, which was relatively unaffected by the coating failure. Although 
this order has a smaller number of lines than the 5225A order, we found that by masking 
out the Ha line when executing the cross-correlation we can get velocities with an error of 
1 — 2 km s^^. This velocity resolution is good enough for our kinematic study (see, e.g., 
Sicilia-Aguilar et al. 2005). While this order is not optimal for radial velocity measurement, 
observations of the Ha profiles allows us to unambiguously identify and distinguish CTTS 
and WTTS stars. 

Our first observations in NGC 2264 were taken in March 2004. These data were obtained 
during the commissioning phase of the instrument, and we ended up observing 155 targets 
in one fiber configuration (see Table [T]) centered north of the Cone nebula. The preliminary 
results derived from this first data set showed interesting results so we decided to observe 
NGC 2264 again. We revisited the cluster selecting targ ets from an expanded list including 



objects from the cataloge of R04 and the full cataloge of iRamirez et al.l (120041 ). 



Even though the spectrograph mirrors were re-coated during the summer of 2004 and 
we could have used the 5225A order, we found that the accuracy of radial velocities obtained 
in the Ha order to be adequate and the extra information in the Ha line profiles to be 
scientifically useful so we continued our observations in this order. The second observing 
run was carried out in December, 2004. During the nights of 2 Dec. and 3 Dec. we observed 
923 targets in four fields (see Tabled]), including 88 targets observed in March which were 
thought to be members based on the March RV data. 

We have found, in observations of regions with strong, spati ally structured nebular 



emission, e.g. the Orion Nebula Cluster (jSicilia-Aguilar et al.ll2004l ). that it is necessary to 



offset the telescope 10" without moving the fiber buttons to obtain a good measurement 
of nebular emission contaminating the stellar spectra. In March we followed this procedure. 
We found that in NGC 2264 this background subtraction is only important near the Ha and 
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[N II] lines; it is unimportant for radial velocity measurements, which depend only upon the 
absorption hne spectra in between the emission features. Therefore in December we opted 
to eliminated the offset observations so that we might spend more time collecting stellar 
spectra. 



2.3. Radial velocities 



To establish the radial velocity zero point for our spectra and to provide a variety 
of well-observed spectra that could serve as templates in the cross-correlation analysis, 
we obtained spectra of 10 re latively bright stars in a field located in Selected Area 57 
(jPickering and Kapteynlll918l ). near the north Galactic Pole. The radial velocities of these 
stars have been monitored for many years using the CfA Digital Speedometers (JLatham 
19921 ) without showing any signs of variation, and their absolute velocities have all been 
measured with a precision better than 0.2 km s~^ following the procedures outlined by 
Stefanik. Latharn. fc Torres! (119991 ). These 10 stars were included in the sample discussed 
by iLatham et al.l ( 199ll ). but two recent technical developments have made it possible to 
improve measurements of their mean velocities. First, several additional observations had 
been obtained over the intervening years, thus doubling the typical time spanned by the 
radial velocities measurements. Second, the spectra were re-analyzed using our latest library 
of synthetic spectral templates calculated by Jon Morse for a grid of iKuruczl (119791 ) model 
atmospheres for a range of effective temperatures, surface gravities and rotational velocities. 
Solar abundance is assumed in all cases. We ran grid s of correlations betw een these tem- 
plates and observed spectra using the rvsao package (JKurtz fc Minklll998l ) running inside 
the IRAF environment. For the final velocity determinations, we adopted the template pa- 
rameters that gave the strongest correlation peak averaged over all the observed spectra (see 
Table [21 and Stefanik et al., in prep.). 

To calibrate the wavelength scale of our spectra, we used Thorium-Argon lamp expo- 
sures for both the NGC 2264 targets and the SA5 7 secondary standards. We adopted the 
observed spectrum of W23870 (JLatham et al.lll99ll ) in SA57 as the template for the velocity 
determinations of all our NGC 2264 targets in March, because this K dwarf gave the best 
results using the cross-correlation technique. However, for the final reduction of the entire 
data set we made use of all the template stars in the SA57 field, in order to obtain the 



-•^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
contract with the National Science Foundation. 
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most accurate radial velocities. For each of the observed NGC 2264 targets we selected the 
"best matching template" which yielded the correlation peak with the highest signal-to-noise 
value. If the RV values from all the other 9 templates were comparable to the radial velocity 
measured with this best matching template, we combined all the 10 measurements of a given 
star and used the median value for the further analysis. 

Using this method of template selection, we obtained radial velocities with errors less 
than ± 1.5 km s~^ for 344 targets, as shown in Fig. [TJ In this figure we also plot the Ha 
emission stars of R04 objective prism survey to guide the eye where the larger concentrations 
of cluster members are expected. 

In order to ensure consistency between the March and December data we observed 88 
stars in both runs. For 31 of these stars yielding the highest correlation peaks the mean of 
the velocity differences is RVMar — RVoec = 0.3 km s~^ with an RMS of 0.8 km s~^. 



2.4. ^^CO observations 

Observations of the ^^CO transition of CO (110.201 GHz) were obtained at the Five 
College Radio Astronomy Observatory (FCRAO) 14m telescope located in New Salem, Mas- 
sachusetts during 2 002 September, as part of a larger survey of nearby cluster forming regions 



(IRidge et al.ll2003f ). 



Using the SEQUOIA 32-element focal-plane array and On-the-Fly (OTF) mapping tech- 
nique eight 15' X 15' submaps were obtained over 3 nights, in order to build up the final 
30' X 60' map. Submaps were obtained by scanning in RA, and an "off" source reference 
scan was obtained after every two rows. The off-position was checked to be free of emission 
by performing a single position-switched observation with an additional 30' offset. Each 
submap was repeated twice to increase the sensitivity. 

The narrow band digital correlator provided a total bandwidth of 25 MHz over 1024 
channels, yielding an effective velocity resolution of 0.07 kms~^. System temperatures at the 
observing frequency were between 200 and 300 K (single sideband). Calibration was via the 
chopper- wheel technique (Kutner & Ulich 1981), yielding spectra with units of T^*. The 
final combined map has an average rms sensitivity of 0.2K/channel. A full discussion of 



the da ta analysis procedure and the integrated intensity image was presented in lRidge et al. 



(|2003|) 
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3. Results 

3.1. Distribution of radial velocities 

The radial velocity distribution of the 344 stars with accurate RV values is shown on 
Fig. [21 The cluster members clearly form a peak centered around V^^^^° = 22 km s^^, but 
the distribution is unusually wide and not symmetrical. We also report the distribution on 
an expanded velocity scale to show how clearly the cluster members form a peak in radial 
velocity space. 

W e plot the RV distribu tion for four selections of stars according to the signal-to-noise 



values (JTonry fc David Il979l ) of the correlation, in Fig. [2l This R value is i nverse ly propor 



tional to the error of velocity measurement. As described in lKurtz fc MinkI (1l998l ) the errors 
consist of a systematic and a statistical component, the systematic could be determined 
for a given instrument by comparing radial velocities measured on large number (several 
hundred) of objects at two epochs. A very good estimate for the random velocity error is 



a = (3/8) * [w/{R + l)] flKurtz fc Minklll998[ ). where w is the FWHM of the correlation peak 



and R is the signal-to-noise ratio. 

This estimate agrees well with our result in NGC 2264 based on 31 stars, observed both 
in March and December with an R value larger than 6 (see §2.3). Also, based on results 
with other instruments, the estimate for the error mentioned above is likely within ±30% 
of the real error value for a Hectochelle measurement. We therefore adopt a velocity error 
of cr = Q/{R + 1) kms"^; the error is larger for stars with rotational velocities greater than 
about 10 km s~^. In principle, this means that slowly-rotating stars with correlation peaks 
of -R = 2 can yield a velocity error of 2 km s~^. 

Of course, in our spectral order of velocity width ~ 8000 km s~^, we can expect to find, 
just due to chance, roughly one 3a peak (with an FWHM of ~ 10 km s^^ or an R value 
of ~ 3) in the cross-correlation function, somewhere within the interval. If we restrict our 
attention to an interval of ~ 100 km s~^ near the median cluster velocity, then we would 
expect to have a probability of about 10% of finding a i? ~ 2 peak of width ~ 10 km s~^. 
The R > 2 selection is included in Fig. [2] because, even with the larger uncertainty in veloc- 
ity, the distribution of these stars follow the overall shape of the histogram fairly well. As a 
reminder that large errors are associated with this group, we plotted this group as light gray 
columns without a border. There is a fifth sample of stars displayed on the left panel of Fig. 
|5] as narrow, white bars: stars with Ha emission detected in their spectra. Comparing this 
sample to the others we can say that the radial velocity ditribution of TTS is the same as for 
the rest of our targets. Also, as the distribution is same independently of the R value of the 
correlation, the observed large dispersion is not due to uncertainties in the RV measurement. 



The RV histogram of Fig. [2] is far from Gaussian. To estimate the velocity dispersion 
we used the definition a^ = X] (^« — x) /{N — 1). Assuming lower and upper limits for 
membership as 13 km s~^and 30 km s~^, the estimated one- dimensional RV dispersion is 
a ~ 3.5 km s~^. This is much larger t han in the ONC, for example, which has a dispersion 



of 2.3 km s~ (1 Jones &: Walked 119881 ). and ~ 1.8 km s~ for the brighter stars (see also 
Sicilia-Aguilar et al. 2005). The skewness of the velocity distribution suggests that the 
cluster is not dynamically relaxed, consistent with our finding of spatially-coherent motions 
(§3.3). 



3.2. Ha emission and binary stars 

Recording spectra at the Ha order has the advantage of not just measuring RV, but 
also identifying T Tauri stars and distinguishing between the CTTS and WTTS. Classical T 
Tauri stars in star forming regions were originally identified as strong emission line objects 
in objective prism spectra. Later other techniques such as X-ray surveys could find similar, 
young stars but with much lo wer optical emission line strengths than in classical TTS. As 



described in iHartmannl (119981 ) WTTS are not accreting from circumstellar disks. In absence 
of this magnetospheric accretion, which would produce emission lines with large velocity 
widths ( > 100 km s~^), solar- type magnetic energy dissipation powers the excess (chromo- 
spheric) emission in WTTS so the observed lines are much narrower and weaker than in the 
case of CTTS. 



We found 277 stars with Ha emission exceeding the nebular compone nt: 152 among 



these were CTTS and 125 were WTTS, using the indicator of accretion from lWhite fc Basri 



( 120031 ). that the Ha full width at 10% intensity exceeded 270 km s~^. A sample of Ha 
profiles plotted on Fig. O We define membership as the RV value should be within 4 a of 
the cluster mean velocity (22 km s~^ ± 4x 3.5 km s~^). Based on this criteria we found 178 
of these emission stars to be radial velocity members (with an R value larger than 2). For 
these stars identification is given in Table [31 A few of the Ha stars listed in Table [3] have 
velocities close to the 4 a RV membership selection criteria, these are noted as hsb. Further 
observations are needed to see if these are spectroscopic binaries. 

Of the 277 obvious cluster member emission stars, 99 (35%) have velocities outside either 
4cr lower or higher than the mean cluster velocity. These stars yielded an R value less than 
2, which suggest a random peak in the cross-correlation happened to be close to the cluster 



-9- 



velocity. Based on their emission spectra these stars are young and therefore members, but 
we do not have enough SNR in the cross-correlation to determine velocities. 

We identified several spectroscopic binary stars. Some of these have resolved double 
peaks in the cross-correlation function, and we measured different velocities in the two sep- 
arate observing runs. There are 6 such stars listed in Table [3] and H] with a three- letter note 
"RDB" (Resolved Double- lined Binary). In case of obvious side-lobes of the cross-correlation 
function and two non-consistent but accurate velocities, we identified additional 3 stars as 
binaries (listed with a note "UDB", as Unresolved Double- lined Binary). There are 6 "rdb" 
and 5 "udb" noted stars in Table [3] and HI for which we have just one RV measurement but 
observed double peaks or side-lobes in the cross-correlation function. 11 stars showed no sign 
of binarity in the cross-correlation, but we got two significantly different velocities. These 
are likely single-lined binaries, and have an "sb" note. 



3.3. Velocity correlations betvi^een stars and gas 

The radial velocity distribution of stars in NGC 2264, as shown in Fig [2l is wide and 
asymmetric, probably due to subcluster groups with different mean velocities. In optical 
images the cluster shows an elongated shape to the north-south direction. The Ha stars 
of R04 also show two denser cores separated mainly to declination, one centered around 
a = 6'^41™10" and 6 = +9°30' and the other around a = 6^40'"50" and 6 = +9°50' (see 
Fig. [1]). To explore this pattern further, in Fig. H] we plot the declination values of stars 
observed at two epochs, as a function of velocity. In order to decrease uncertainties due to 
measurement errors, and natural distribution in velocity space within a possible subcluster 
group, we calculated mean declination values for 2 km s~^ wide bins in RV (Fig. H]). The 
shaded areas show the declination ranges of the two largest (suggested) subclusters based 
on the surface density of Ha emission stars in the catalog of R04. The southern-most 
declination limit of the shading was set by the Hectochelle field observed in March. The 
increasing velocity with increasing declination is obvious. 

There is not a clear trend in RA, but this is consistent with the optical image of the 
cluster which shows a North- South elongation with only a very slight tilt in RA. But the 
spatial pattern of stellar velocities becomes clear if we compare them to the pattern of molec- 
ular gas. Using ^^CO observations we tried to match the motion of the gaseous component 
with the motion of the stellar component in NGC 2264. Fig. [5] displays channel maps as 
eight, 1 km s~^ wide panels centered between Vlsr = 3-10 km s~^. Stars from the peak of 
Fig. [2] with the same 1 km s~^ binning are superimposed on these maps; the mean velocity 
of these groups are VheUoc = 17-24 km s~^. The correlation between the position of dense 
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CO cores and position of stars is very convincing. 



On the upper-right panels of Fig. \5\ there is a dense ^^CO core at a = 6''41™00^ and 
6 = +9°36'. On the same panels, many stars occupy the same region (B) in the given RV 
bin. As we go to higher RV values, stars disappear from this area, as well as the gas, but 
two other dense ^^CO regions appear. One is to the north, with an east-west elongation and 
a complicated structure, around the central point of a = 6'*40'"40^ and 6 = +9°52'. Most of 
the stars fall in this region (A) on the lower panels of Fig. [5l especially on the lower-right ones. 

On the lower-middle panels of Fig. [5l indicated by "C", there is another dense ^^CO 
core at a = 6'^41'"10'* and 6 = +9°29' at these higher velocity values. This surrounds the 
Cone-nebula. On these panels, there are several stars very close to this region, well separated 
in space from the larger number of stars to the north but lieing in the same radial velocity 
range. These stars are responsible for the turn-down of the north-south gradient on Fig. |H 
This can be a hint of a third subcluster group around the Cone, having a mean heliocentric 
velocity of ~ 22 km s~^. 

The distribution of R04 Ha stars also show correlations with the densest cores of 
molecular gas. At Vlsr = 9 km s~^ there is a small, 3' diameter ring of stars around 
a = 6''41™10^ and 6 = +9°29', which surrounds a strong concentration of molecular gas 
best seen on that panel. The ring-like distribution of RV measured member stars on the 
Vlsr = 4 km s~^ panel could also be a sign of some high energy process which played role 
in the formation of the cluster. 

Similar agreements in the structure of stellar and gaseous components can be seen on 
Fig. O which is a projection in RA, displaying the velocities along the horizontal axis and 
DEC values along the vertical axis. At Vlsr = 7 km s~^ and 6 = +9°53' there is a"hole" 
in the gas, likely a projection of a bubble, surrounded by stars in a ring which has higher 
stellar density than the neighboring regions. This structure likely has a connection with the 
star HD 47839 at a = 6^40"^58" and S = +9°53'44" 



4. Discussion 

While other star-forming regions have been shown to have significant substructure in 
the young star population (e.g., Lada & Lada 1995; Gomez et al. 1993), NGC 2264 is the 
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first region for which we have sufficiently accurate measurements to explore substructure in 
the stellar radial velocities. The spatial distribution of stellar radial velocities and its corre- 
lation with the velocity structure of the molecular gas shows that NGC 2264 is composed of 
relatively distinct components (Figures EED- Intrinsic errors due to various types of activ- 
ity (starspots, accretion; Guenther et al. 2000) are not large enough to obscure the strong 
correlation of the stellar velocities with the gas (Figures [5]|6|) . Though individual stars could 
have velocity offsets on the order of few krn s~^ due to fast rotation and starspots/accretion 
(JNeuhaeuser et al.lll998l : IStrassmeier et al.ll2005l ). that is very extreme and the large number 
of stars in our sample overrides the significance of such outliers. 

The presence of such substructure in NGC 2264 is not surprising given its size and youth. 
Assuming a local, within-group velocity dispersion (one-dimensional) of order 2 km s~^, and 
a spatial extent of ± 4 pc (~ ±17 arcmin) or so from the central region, the resulting crossing 
time of ~ 2 Myr is comparable to the estimated age of the cluster (see Introduction). Even 
adopting a velocity dispersion of 3.5 km s~^ results in a crossing t ime of ~ 1 Myr , close 



to the youngest stellar age estimates. Numerical simulations show (jScally fc Clarkd l2002l ) 
that substructure in a stellar cluster can remain detectable even after several crossing times, 
depending upon the number of stars involved and their initial subcluster densities. However, 
it is unlikely that the region is considerably older than one or two crossing times, because the 
gas also shows the same substructure as the stars. The cluster gas will shock and dissipate 
considerable amounts of energy on the ffist crossing, and thus the initial substructure should 
be damped out much more rapidly than for the stellar population. 



Crutcher fc Hartkopil (119781 ) presented observational evidence for two separate gaseous 
cores inside NGC 2264 which have a mean velocity difference of 2 km s~^. These two cores 
correspond in position with the centers of our southern subclusters B and C ( C is in the 
direct vicinity of the Cone nebula), see lower left panel on Fig. O The velocity gradient we 
est imate based on our ob s ervati o ns agrees in direct i on an d in its amount with that determined 
by lCrutcher fc Hartkopil (119781 ). iLang &: Willsonl (Il980l ) verified the same gradient based on 
NH3 measurements and also explained the motion of the two cores as rotation, with a period 
of 4 Myrs if th e rotational axis is perpendicular to the line of sight. However, as the ^^CO 
observations of iRidge et al.l (120031 ) show (Fig. E]), the velocity field is more complicated than 
simply that of pure rotation. 

Large velocity dispersions sometimes arise in molecular gas as the result of powerful 
bipolar outflows or other stellar energy input. However, the systematic motions seen in the 
gas require too much energy to explain in this way. Instead, we interpret the velocities as 
mostly gravitationally-generated, which would very simply explain why the stars spatially 
co-located with the gas show the same velocities. 
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Burkert &: Hartmannl (120041 ) recently explored a simple model of the formation of star- 
forming filaments. Their calculations assumed that molecular clouds were formed schemati- 
cally from colliding supersonic flows which formed (approximately) sheet-like configurations. 
If the sheet is elongated, the result is gravitational collapse first along the shortest dimen- 
sion, forming a filamentary structure. As shown in the left panel of Fig. [7] (reproduced 
from Figure 10 in Burkert & Hartmann 2004), irregularities initially present in the sheet 
grow to form non-linear density enhancements by the time the configuration collapses to 
a filament. The right panel of Fig [7] shows the integrated surface density as a function of 
the radial velocity in the x-direction, showing that multiple dense structures with coherent 
velocities are present. Often individual "clumps" appear in nearby pairs with distinct ve- 
locities. This pattern is the result of the overall collapse of the elongated sheet; material 
originating on opposite sides of the sheet fall in with opposing velocities. The high- velocity 
infall at each end of the filament is due to the highest gravitational acceleration at filament 



ends (JBurkert fc Hartmannl 120041 ). coupled with some initial rotation of the sheet. 



We suggest that the Burkert & Hartmann model qualitatively explains the velocity 
structure seen in the NGC 2264 region. The pattern of first higher velocities, then lower 
velocities, then somewhat higher again proceeding southward along the region (Figure 7) can 
be explained by gravit at ionally- driven infall from opposite sides of a cloud falling in to form 
a more filamentary structure. The highest velocity structure seen at the northern end may 
be similar to the acceleration seen in the end of the filament in the right panel of Fig. [7l The 
velocity uni t s in F igure 7 are (GMR/nY^'^; if we adopt a total gas mass of 4OOOM0 from 



Ridge et al.l (120031 ) and a filament half-length of 4 pc, these units correspond to ~ 2 km s^^. 
Clearly, gravitational collapse can account for the magnitude of the motions seen in NGC 
2264, especially considering geometric uncertainties and the simplicity of the model. 



5. Summary 

We have carried out a spectroscopic study of stars in NGC 2264, utilizing the Hec- 
tochelle multiobject spectrograph on the 6.5 m MMT. Obtaining 1078 spectra of 990 stars 
we classified 471 stars as cluster members based on their radial velocity and/or Ha emis- 
sion. We determined radial velocities with an accuracy of < 1.5 km s~^ for 344 members. 
The radial velocity distribution of these stars is wide, non-gaussian, and the dispersion is 
3.5 km s~^. We compared the stellar radial velocities to the velocity of the molecular gas in 
the cluster, traced by ^^CO, and found a strong correlation in the spatial distribution of the 
two components. At least three subgroups of stars can be identified with distinct spatial and 
velocity coherence. We interpret this substructure as the result of gravitational collapse of 
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initial clumps of star-forming gas from a more extended structure to a roughly filamentary 
distribution. Using the results presented here, more advanced numerical simulations could 
help elucidate the initial conditions that produced NGC 2264. 

This work was supported in part by NASA grant NAG5-13210. 
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Fig. 1. — Our targets in the field of NGC 2264 (open diamonds), superimposed on a DSS 
POSSl red plate image of the cluster. Ha emission stars (filled circles) from R04 are also 
shown to guide the eye where the concentration of young stars expected to be larger. Most 
of these emission stars were observed by us, but as can be seen our targtes were more evenly 
distributed, covering the outer regions of the cluster as well. 
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Fig. 2. — Radial velocity distribution for 344 of the observed stars which were giving con- 
sistent velocities independent of the RV template. The distribution is unusually wide and 
non-gaussian. The insert is an expanded velocity range to show how clearly the cluster mem- 
bers form a peak in the velocity space. There are four selections of stars plotted according 
to the signal-to-noise (R) values of the correlation, showing that the large dispersion is real 
and not coming from the uncertainties of RV measurements (see text for details). The white 
impulses show a selection of stars with Ha emission presented in their spectra. The similar- 
ity between the distribution of these and non emitting stars suggest that even a selection is 
biased towards TTS it can be representative of the entire cluster in kinematical studies. 
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Fig. 3. — Examples of Ha profiles for stars with detected emission. Wavelength is converted 
to RV, and the 2MASS identification number is shown in each stamp. The emission of the 
background is not subtracted, as we did not have separated offset sky exposures for all of 
our fields. For all profiles see the electronic format of the paper or contact the authors. 
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Fig. 4. — A north-south velocity gradient observed in NGC 2264 is clearly shown, by plotting 
the mean declination values in 2 km s~^ wide RV bins against the mean radial velocity of the 
bins. Only stars observed at two epochs (and giving consistent velocities) used for this plot, 
in order to decrease scatter. The error bars represent the RMS of DEC/RV values of stars 
in a given velocity bin. The small DEC error bar for the RV group of 24-26 km s~^ means 
that almost all the stars in that velocity bin are at the northern-most part of the cluster. 
The shaded areas show the declination ranges for the condensations of Ha stars in R04. 
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Fig. 5. — We find strong spatial correlation between the molecular gas and stars in velocity 
space, by plotting radial velocity members of our sample (filled boxes) superimposed on ^^CO 
channel maps (both averaged in 1 km s~^ wide bins in RV). Positions of Ha emission stars 
from R04 are also shown as dots to draw an outline of the cluster (see upper left panel). 
Three supposed subcluster structures (A,B,C) are also shown in the RV channels in which 
those are prominent. 
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Fig. 6. — ^'^CO spatial-velocity plot of NGC 2264 projected along a shows an increased ve- 
locity dispersion at the northern end of the filament, a north-south gradient similar observed 
in the stellar component also plotted on Fig. HI The molecular gas is shown as a shaded 
area, the dots are member stars with measured velocities. The size of the dots correspond to 
the accuracy of RV measurement (bigger dot means higher accuracy) . Note the bubble-like 
feature at V=7 km s~^, 6 = +9°53', surrounded by stellar members. This structure likely 
connected to the bright (V=4.6) star HD 47839 at a = 6''40™58" and 6 = +9°53'44" 



-22- 



i~ 0.0 




-0.5 



50.0 




,''''■''''' 


41.7 




* 




5 


_ 


33.3 




II 




> 


^^- 


■ 


Z5.0 


■°-'*^sMf 


h 




^^"Vl 


n 


16.7 




^H 1 




-5 


1 


B.3 




» 


0.0 


-10 


' 



IB. 7 



l5.a 



2.9 



11.5 



-1.0 



-QJi 



0,0 



0.5 



Fig. 7. — Simulation of gaseous sheet collapse to a filament (see Burkert & Hartmann 2004). 
The left panel shows the spatial distribution of the filament, which exhibits a number of 
condensations which resulted from slight substructure in the initial conditions, subsequently 
magnified by gravitational collapse from an initially elliptical configuration. The right panel 
shows the projected emission along one axis (assumed optically-thin) as a function of radial 
velocity. The various condensations show offset velocities which are the result of gravitational 
collapse from opposite sides of the filament. The scaling to the mass and radius of the NGC 
2264 gaseous filament is such that each unit of velocity corresponds to 1.6 km s~^. The 
existence of condensations with differing velocities at the ~ 2 km s~^level, plus the higher- 
velocity tails at the end of the filament - the result of gravitational attraction on at "focal 
points" (see Burkert & Hartmann 2004) where especially large density concentrations develop 
- is suggestive of the structure seen in both the gas and stars of NGC 2264 (Fig. [6]) (see 
text). 



Table 1. Summary of spectroscopic observations 



Date 


field center 




exposures 


off'set sky 


number of targets 


2004 Mar 10 


Qh 4om ^^s 


+ 09° 44' 


46" 


4 X 15 min 


1 X 15 min 


155 


2005 Dec 02 


Qh 4om 4gs 


+ 09° 38' 


13" 


3 X 15 min 


— 


234 


2005 Dec 02 


6*^ 40"^ 35= 


+ 09° 36' 


56" 


3 X 15 min 


— 


232 


2005 Dec 02 


gh 4Qm 42S 


+ 09° 45' 


36" 


3 X 15 min 


— 


231 


2005 Dec 03 


Qh 4om 42S 


+ 09° 32' 


13" 


3 X 15 min 


— 


226 



to 

00 



Table 2. Observed radial velocity standard stars in a sub-field of Selected Area 57 



Name 


RA 


DEC 


Vrad 


Vrot 


Teff 


log{g) 


W33245 


13:03:50.9 


+30:11:38 


-20.95 


0.0 


5000 


3.5 


W22942 


13:04:25.2 


+30:25:44 


-16.32 


3.3 


4500 


2.5 


W24128 


13:04:45.0 


+29:51:22 


2.20 


3.2 


5500 


4.5 


W23082 


13:04:58.4 


+30:10:39 


-21.68 


3.5 


6250 


4.5 


W23108 


13:05:18.4 


+30:10:46 


-16.43 


7.0 


4250 


2.0 


W24226 


13:05:23.5 


+29:43:02 


11.69 


4.1 


6500 


4.5 


W23131 


13:05:29.3 


+30:11:12 


-5.96 


5.4 


5500 


3.5 


W23961 


13:06:00.6 


+30:01:22 


-6.98 


4.2 


6000 


4.5 


W23870 


13:06:45.8 


+29:46:30 


-14.85 


5.3 


4750 


4.0 


W23833 


13:06:58.6 


+29:57:49 


-4.90 


0.5 


5250 


4.0 



Note. — Table Vrad — heliocentric velocities on the native CfA velocity 

system. To convert to an absolute velocity system based on extensive |.^ 

observations of minor planets, 0.139 km s~ ^ should bo added to the native nf^ 
CfA velocities (iStefanik. Latham, fc Torres[[1999i) . i 



Table 3. Members 



2MASSJd 


REBU REIP LAMM J 


(J^H) ( 


H -K) 


Ha 


RVoec/RVMar 


Rdcc/ Rlvlar 


notes 


0639037+094023 


325 


- 12.65 


0.74 


0.35 


c 


20.46/- 


7.6/- 


H 


0639156+093912 


- 


- 13.26 


0.44 


0.11 


- 


17.16/- 


4.0/- 


H 


0639157+092929 


- 


- 11.67 


0.69 


0.24 


- 


33.80/- 


12.3/- 


H 


0639170+095815 


- 


- 12.43 


0.36 


0.17 


- 


14.06/- 


6.2/- 


H 


0639226+091703 


- 


- 11.13 


0.19 


0.14 


- 


28.28/- 


2.2/- 


H 


0639226+095430 


330 


- 14.25 


0.61 


0.35 


w 


16.51/- 


2.1/- 


H 


0639229+093616 


- 


- 12.62 


0.49 


0.14 


- 


14.49/- 


11.3/- 


H 


0639246+092233 


- 


- 12.79 


0.41 


0.12 


- 


18.30/- 


9.7/- 


H 


0639262+094110 


R453 


- 12.71 


0.58 


0.23 


w 


16.17/- 


6.2/- 


H 


0639266+093037 


- 


- 15.50 


0.66 


0.25 


- 


29.76/- 


2.5/- 


H 


0639275+092309 


- 


- 11.03 


0.91 


0.31 


- 


34.54/- 


11.2/- 


H 


0639278+092839 


- 


- 10.81 


0.54 


0.20 


- 


13.17/- 


11.6/- 


H 


0639293+092058 


- 


- 12.20 


0.41 


0.12 


- 


25.38/- 


11.2/- 


H 


0639303+093337 


- 


- 12.93 


0.69 


0.14 


w 


21.27/- 


3.6/- 


H 


0639329+095629 


R691 


- 13.16 


0.46 


0.13 


- 


23.95/23.68 


10.7/9.4 


H 


0639333+095201 


R705 


- 12.63 


0.73 


0.24 


c 


19.32/- 


4.1/- 


H 


0639335+095140 


R712 


- 13.71 


0.30 


0.30 


- 


25.61/25.83 


3.1/3.0 


H 


0639339+094920 


R729 


- 13.04 


0.77 


0.17 


w 


19.92/- 


5.3/- 


H 


0639344+095451 


- 


- 12.99 


0.78 


0.28 


c 


19.25/- 


6.4/- 


H 


0639347+094654 


- 


- 14.45 


0.59 


0.26 


- 


13.30/- 


2.8/- 


H 


0639353+093232 


- 


- 13.21 


0.44 


0.14 


- 


28.06/- 


2.9/- 


vm 


0639360+092426 


- 


- 11.85 


0.67 


0.15 


- 


13.90/- 


11.6/- 


H 


0639388+095151 


R893 


- 12.47 


0.40 


0.23 


- 


13.61/- 


3.0/- 


H 


0639393+094521 


R908 


- 12.90 


0.86 


0.50 


c 


19.26/- 


7.8/- 


H 


0639393+095559 


R905 


- 12.79 


0.65 


0.24 


w 


18.32/- 


4.2/- 


- 


0639414+094619 


R985 


- 12.57 


0.65 


0.22 


c 


18.95/- 


7.5/- 


H 



to 



Note. — Hectochelle targets in NGC 2264 found to be members based on measured RV value or detected Ha emission. The criteria 
of becing RV member is to have at least one velocity measurement within 4a- of the cluster mean velocity: 8 km s~^ < Vhelio < 36 
km s^-*^ . (Stars listed only with R value of cross correlation larger than 2.) For the entire table see the electronic version of the Jurnal, 
or contact the authors. 

2MASS id — identification number from 2MASS catalog (truncated RA and DEC coordiantes as: HHMMSSS+DDMMSS); REBU — 
id numbers from RebuU et al 2002, AJ123, 1528; REIP — id from Reipurth et al 2004 AJ127, 1117; LAMM — id from Lamm et al 



2004 A&A 417, 557; J — 2MASS J magnitude and; (J-H) — 2MASS (J - H) color index; (J-H) — 2MASS (H ~ K) color index; 
Ha — Ha emission detected: c - CTTS; w - WTTS; x - V642 Mon; r - reduction error in 2004 Dec data set, categorized as WTTS in 
2004 Mar data; 

VR£)ec/VR,jv/ar — heliocentric radial velocity measured 2004 Dec / 2004 Mar; ^uec/^Mar — R- value of cross correlation (see text for 
details) 2004 Dec / 2004 Mar; 

notes — MEMBERSHIP: vm - at least one RV measurement suggest membership, but R value is low; hm — just Ha emission suggest 
membership, none of the RV values (if there are two); ACCURACY: H - most of the templates gave the same velocity, R value of 
cross-correlation is high, accurate RV value (typical error < 1.5 km s^^) for one RV measurement; BINARITY: RDB — resolved, double 
lined binary, with RV at 2 epochs (separated double peak in correlation function); UDB — unresolved double lin. binary, with RV at 2 
epochs (blended/side lobed peak, but clear asymmetry in correlation function); rdb — same as for RDB, but only 1 RV measurement; udb 
- same as for UDB, but only 1 VRV measurement; sb - likely single lined binary, RV at two epochs and at least one suggest membership, 
some also have Ha emission; hsb - member based on Ha emission detected, but correlation either had just acceptable R value (barely 
larger than 2) or radial velocity is close to the edge of RV distribution. 



to 



Table 4. Non-members 



2MASSJd 


REBU REIP 


LAMM J 


(J^H) { 


H -K) Hi 


0, RVoec/RVMar 


Roec/RMar 


notes 


0639152+093937 


R3 


- 10.90 


0.76 


0.27 


45.53/- 


5.5/- 


H 


0639161+092105 


- 


9.95 


1.00 


0.36 


85.84/- 


11.0/- 


- 


0639163+091953 


- 


- 12.61 


0.29 


0.15 


65.74/- 


2.8/- 


- 


0639174+092137 


- 


- 13.61 


0.32 


0.18 


44.93/- 


2.3/- 


H 


0639183+093336 


- 


- 14.88 


0.71 


0.13 


42.58/- 


2.7/- 


H 


0639200+093036 


- 


- 11.32 


0.70 


0.24 


49.10/- 


5.2/- 


H 


0639229+093523 


- 


- 12.41 


0.45 


0.13 


-12.36/- 


9.0/- 


H 


0639239+094201 


R353 


- 10.84 


1.09 


0.44 


81.66/- 


11.3/- 


- 


0639240+092550 


- 


- 11.77 


0.67 


0.19 


44.22/- 


5.4/- 


H 


0639241+091805 


- 


- 13.48 


0.37 


0.15 


-19.37/- 


5.0/- 


H 


0639243+095619 


- 


- 10.68 


0.74 


0.28 


-12.91/- 


11.6/- 


H 


0639245+095040 


- 


- 10.63 


0.26 


0.07 


85.41/- 


10.4/- 


- 


0639246+091721 


- 


- 12.84 


0.81 


0.34 


69.63/- 


6.2/- 


- 


0639250+095339 


- 


- 13.06 


0.51 


0.14 


78.26/- 


3.9/- 


rdb 


0639275+095100 


- 


- 11.67 


0.79 


0.28 


87.58/- 


12.1/- 


- 


0639283+095358 


- 


- 10.06 


0.82 


0.31 


0.78/- 


7.5/- 


H 


0639287+092640 


- 


- 12.40 


0.28 


0.16 


-2.97/- 


2.5/- 


H 


0639303+092300 


- 


- 12.05 


0.47 


0.10 


109.22/- 


10.9/- 


- 


0639304+094924 


- 


- 10.03 


0.77 


0.25 


-9.34/- 


12.4/- 


H 


0639311+095330 


- 


- 10.32 


0.28 


0.08 


4.17/- 


9.8/- 


H 


0639333+093011 


- 


- 13.56 


0.47 


0.17 


64.57/- 


4.1/- 


- 


0639334+094728 


R707 


- 13.36 


0.40 


0.20 


37.26/- 


2.0/- 


H 


0639344+092338 


- 


- 12.13 


0.38 


0.10 


99.31/- 


4.8/- 


- 


0639346+092440 


- 


- 12.72 


0.35 


0.07 


52.81/- 


5.5/- 


- 


0639373+092637 


- 


- 11.47 


0.63 


0.17 


39.21/- 


10.2/- 


H 


0639386+094044 


- 


- 12.74 


1.27 


0.57 


52.73/- 


2.0/- 


- 


0639396+094544 


R916 


- 11.98 


0.62 


0.19 


-/43.05 


-/7.2 


- 



to 



Note. — Hcctochellc targets in NGC 2264 found to be non members based on RV measurement. (Stars listed only with R value of 
eross eorrelation larger than 2.) For the entire table see the electronic version of the Jurnal, or contact the authors. 

2MASS id — identification number from 2MASS catalog (truncated RA and DEC coordiantes as: HHMMSSS+DDMMSS); REBU — 
id numbers from RebuU et al 2002, AJ123, 1528; REIP — id from Reipurth et al 2004 AJ127, 1117; LAMM — id from Lamm et al 
2004 A&A 417, 557; J — 2MASS J magnitude and; (J-H) — 2MASS (J - H) color index; (J-H) — 2MASS {H - K) color index; 



Ha — Ha emission detected: c - CTTS; w - WTTS; x - V642 Mon; r - reduction error in 2004 Dec data set, categorized as WTTS in 

2004 Mar data; 

VR£)ec/VRMar — heliocentric radial velocity measured 2004 Dec / 2004 Mar; R^jec/R-Mar — R value of cross correlation (see text for 

details) 2004 Dec / 2004 Mar; 

ACCURACY: H - most of the templates gave the same velocity, R value of cross-correlation is high, accurate RV value (typical error 

< 1.5 km s~ ) for one RV measurement; BINARITY: RDB - resolved, double lined binary, with RV at 2 epochs (separated double peak 

in correlation function); UDB - unresolved double lin. binary, with RV at 2 epochs (blended/side lobed peak, but clear asymmetry in 

correlation function); rdb — same as for RDB, but only 1 RV measurement; udb - same as for UDB, but only 1 VRV measurement; sb 

- likely single lined binary, RV at two epochs and at least one suggest membership, some also have Ha emission; hsb - member based 

on Ha emission detected, but have unusable RV measurement due to low SNR; only three of these have trustable velocities with R > 2.5 

(see text). 



to 
00 



